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The effect of quenched disorder on the low-energy properties of various antiferromagnetic spin 
ladder models is studied by a numerical strong disorder renormalization group method and by den- 
sity matrix renormalization. For strong enough disorder the originally gapped phases with finite 
topological or dimer order become gapless. In these quantum Griffiths phases the scaling of the 
energy, as well as the singularities in the dynamical quantities are characterized by a finite dynam- 
ical exponent, z, which varies with the strength of disorder. At the phase boundaries, separating 
topologically distinct Griffiths phases the singular behavior of the disordered ladders is generally 
controlled by an infinite randomness fixed point. 
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I. INTRODUCTION 

Low dimensional quantum spin systems, chains and 
ladders are fascinating objects, which are the subject 
of intensive experimental and theoretical research. The 
main source of this activity is due to the observation that 
quantum fluctuations could result in qualitatively differ- 
ent low-energy behavior in these interacting many-body 
systems. It was first Haldaned who conjectured that an- 
tiferromagnetic (AF) spin chains with integer spin have 
a gap in the energy spectrum (Haldane phase), whereas 
the spectrum of chains with half-integer spins is gapless. 
By now a large amount of experimental and theoretical 
evidence have been collected in favor of the Haldane con- 
jecture. ItJias been realized by Affleck, Kennedy, Lieb 
and Tasaki3 (AKLT) that the ground state structure of 
the Haldane phase for S = 1 is closely related to that of 
the valence-bond solid model, where the ground state is 
built up from nearest-neighbor valence bonds. The hid- 
den, topological order in the chain is measured by the 
non-local string order parameter^: 
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and 



where S" is a spin-1 operator at site i, a = x,y,z 
(. . .} denotes the ground state expectation value. 

Another source of activity in the field of low- 
dimensional quantum spin systems is due to the discov- 
ery of spin ladder materials!! It has been realized that 
spin ladders with even number of legs have a gapped 
spectrum, whereas the spectrum of odd-leg ladders is 



gaplessi. For two-leg ladders, which are analogous ob- 
jects to S = 1 spin chains, the ground state structure 
can be related to nearest-neighbor valence bonds and a 
topological hidden order parameter, similar to that in 
Eq.(l.l) can be defmedEl. 



More recently, ladder models with competing interac- 
tions, such as with staggered dimerizationO and with rung 
and diagonal couplings^, have been introduced and stud- 
ied. In these models, depending on the relative strength 
of the couplings, there are several gapped phases with 
different topological order, which are separated by first- 
or second-order phase transition lines. 

Disorder turns out to play a crucial role in some ex- 
periments on low dimensional magnets. For instance, 
the NMP-TCNQ compoundB can be well described by 
5 = 1/2 spin chains with random AF couplings. More 
recently, non magnetic substitutions in low dimensional 
oxydes such as CuGepaErii 2 ] (being a spin-Peierls com- 
pound), PbNigyj^Og" (being a Haldane gap compound) 
or Y2BaNi05L 4 fO (being a Haldane gap compound) have 
been the subject of intense investigations. The essential 
feature of these compounds is the appearance of anti- 
ferromagnetism at low temperature which can be well 
described by the effective low energy models introduced 
in Refs. 19 Sr(Cui_ x Zn x )2 O3 is a realization of the 



two-leg ladder, and can be doped by Zn, a non magnetic 
ioncj. The specific heat and spin susceptibility exper- 
iments indicate that the doped system is gapless even 
with low doping concentrations. We note that the exper- 
imentally found phase diagram of this compound, as well 
as other quantities, such as staggered susceptibility h^vc 
been obtained by quantum Monte Carlo simulationsEfJ. 
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Theoretically, spin chains in the presence of strong dis- 
order can be conveniently studied by a real-space renor- 
malization group (RG) method introduced by Ma, Das- 
gupta and HuCJ (MDH) . In this method strong bonds in 
the system are successively eliminated and other bonds 
are replaced by weaker ones through a second ordet-per- 
turbation calculation. As realized later by Fishein3 for 
the random spin-1/2 chain and for the related model 
of random transverse-field Ising spin chaii£3 the prob- 
ability distribution of the couplings under renormaliza- 
tion becomes broader and broader without limit and 
therefore the system scales into an infinite randomness 
fixed point (IRFP), where the MDH renormalization be- 
comes asymptotically exact. Fisher has also succeeded 
to solve the fixed-point RG equations in analytical form 
and to show, that for any type of (non-extremely sin- 
gular) initial disorder thjej-system scales into the same 
IRFP. Later numericalcaEl and analytical!!] work has 
confirmed Fisher's results. 

Generalization of the MDH approach for AF chains 
with larger values of the spin is not straightforward, since 
for not too strong initial disorder the generated new cou- 
plings could exceed the value of the already decimated 
ones. To hancUe, this problem for the S = 1 chain Hy- 
man and Yaiigp and independently Monthus, Golinelli 
and Jolicaeuit!3 have introduced an effective model with 
spin-1 and spin-1/2 degrees of freedom and with random 
AF and ferromagnetic (F) couplings. From an analysis 
of the RG equations they arrived to the conclusion that 
the IRFP of the model will be attractive if the original 
distribution parametrized by the power-law form 
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is strongly random, i.e. if < D^ 1 < D^ 1 . For weaker 
initial disorder the system is still gapless, which is called 
the gapless Haldane phase. 

Theoretical work about disordered spin ladders is 
mainly concentrated on the weak disorder limit. Results 
in this direction are obtained in the weak interchain cou- 
pling limit via the bosonization appcaacho and by the 
random mass Dirac fermion methodt-X In particular a 
remarkable stability of the phases of the pure system 
against disorder with XY symmetry has been observeaH. 

In the experimental situation, however, as described 
before the effect of disorder is usually strong and we 
are going to consider this limiting case in this paper. 
Our aim is to provide a general theoretical background 
for strongly disordered spin ladders by studying in de- 
tail several models (conventional ladder, dimerized lad- 
der, zig-zag ladder, and the full ladder with rung and 
diagonal couplings), which could have experimental rele- 
vance. Since often a small change in the couplings or in 
the strength of disorder could cause large differences in 
the low-energy singular properties of the models, there- 
fore we have studied the phase diagrams in the space of 
several parameters. As a method of calculation we used a 
numerical implementation of the MDH approach, which 



could treat the combined effect of disorder, frustration, 
correlations and quantum fluctuations, while some prob- 
lems are also studied by density matrix renormalization 
(DMRG). In particular we have investigated the stability 
of the different topologically ordered phases and studied 
the region of attraction of the IRFP. 

The structure of the paper is the following. In Sec. 
2 we define different spin ladder models and present 
their phase diagram for non-random couplings. A short 
overview about the MDH RG method and its application 
to random spin chains are given in Sec. 3. Our results 
about random spin ladders are presented in Sec. 4 and 
discussed in Sec. 5. 



II. THE MODELS AND THEIR PHASE 
DIAGRAM FOR NON-RANDOM COUPLINGS 

We start with two spin-1/2 Heisenberg chains, labeled 
by t = 1, 2 and described by the Hamiltonians 
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where S/ iT is a spin-1/2 operator at site / and on chain r 
and Ji_ T > 0. For non-random spin chains dimerization 
can be introduced as 
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with n(r) = 0, 1, whereas for random dimerized couplings 
the even and odd bonds are taken from different distri- 
bution. The pure chain without dimerization (7 = 0) has 
a gapless spectrum, and spin-spin correlations decay as 
a power for large distance, which is called as quasi-long- 
range-order (QLRO). Introducing dimerization for 7 > 
a gap opens in the spectrumEj, which is accompanied by 
non- vanishing dimer order, 0^ im 7^ 0. This is measured 
as t he difference between the string order parameters in 
Eq(l.l) calculated with spin-1/2 moments at even (e) 
and odd (o) sites: 
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In the following we generally consider non-dimcrizcd 
chains, otherwise it is explicitly mentioned. 
Now we introduce the interchain interaction 
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which describes the usual rung coupling between the lad- 
ders (see Fig. |l|a). The conventional ladder model is de- 
scribed by the Hamiltonian: H = H 1 + H2 + Hr. In 
the pure model by switching on the AF rung couplings, 
Jp = J R > 0, a Haldane- type gap opens above the 
ground state and the system has a non-vanishiag even 
string topological order, which is measured bytffiJ: 
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For strong AF rung couplings every spin-pair on the 
same rung form a singlet, therefore this phase is called 
the rung singlet (RUS) phase. 
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FIG. 1. Spin- ladder models used in the paper. The con- 
ventional two-leg ladder (a) and staggered dimerization in the 
chain couplings (b). The zig-zag ladder (c) and its represen- 
tation as a chain with first and second neighbor couplings (d). 
The full ladder with rung and diagonal couplings (e). 

Dimerization of the chain couplings could occur in two 
different ways. For parallel dimerization, when equal 
bonds in the two chains are on the same position, i.e. 
in Eq.(2.2) n(l) = n(2), the combined effect of rung 



coupling and dimerization will always result in a gapped 
phase. In the other possible case of staggered dimeriza- 
tion, i.e. with n(l) = — n(2) (see Fig. |l].b), the two chains 
have an opposite dimer order, which competes with the 
rung coupling. As a result the phase diagram of the sys- 
tem (see Fig. ||) consists of two gapped phases, which 
are separated by a gapless transition line, starting in the 
pure, decoupled chains limitQ. 




FIG. 2. Schematic phase diagram of the two-leg AF lad- 
der with staggered dimerization (see Fig. |l|.b for the definition 
of the couplings.) At the phase boundary between the rung 
singlet and dimer phases the gap vanishes. 

Next, we extend our model by diagonal interchain cou- 
plings, given by the Hamiltonian term: 



Hz = 2_^J\ Sj.aSj+iji 
i=i 
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The complete Hamiltonian, H = Hi + H 2 + Hr + Hz, 
describes a zig-zag ladder (see Fig. |l].c) or can be consid- 
ered as a spin chain with nearest neighbor ( J/ 1 , J^) and 
next-nearest neighbor (J/) couplings (Fig. |l|.d). The pure 



model with J, R 



Ji and Ji — J2 has two phases: a 



gapless phase for J2/J1 < -24 is separated from a gapped 
phase by a quantum phase transition point. 

Finally, we extend our model by two types of diagonal 
couplings, which are represented by the Hamiltonian: 
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It is known that the pure AF diagonal ladder described by 
the Hamiltonian, H =J1 1 + H 2 + H D with jf = J D > 
has a gapped spectrumEl. Its ground state is of the AKLT- 
type and has a no n- va nishing odd string orderEj, defined 
in analogy to Eq.(2.5) 
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In the full ladder there are both rung and diagonal 
couplings (see Fig. |].e) and it is described by the Hamil- 
tonian H = Hi + H 2 + Hr + Hjj ■ For non-random AF 
couplings there is a competition between rung and diago- 
nal couplings, so that the ground state phase diagram of 
the system consists of two topologically distinct gapped 
phases (see Fig. ||). The .phase transition between the 
two phases is of first orderB. 
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FIG. 3. Schematic phase diagram of the full AF ladder 
with homogeneous rung and diagonal couplings. The transi- 
tion between the two topologically distinct gapped phases is 
of first order, except of the limit Jr = Jd = 0. 

The main subject of our paper is to investigate how 
the phase diagrams of the pure ladder models, in partic- 
ular in Figs. |2| and [| are modified due to the presence of 
quenched disorder. 



III. THE MDH RENORMALIZATION: RESULTS 
FOR SPIN CHAINS 

In the MDH renormalization group method for ran- 
dom spin-1/2 chains the random AF bonds are arranged 
in descending order according to their strength, and the 
strongest bond, say J23, connecting sites 2 and 3 sets the 
energy scale in the problem, Q = J23. We denote by 1 
the nearest neighbor site to 2 with a connecting bond J12 
and similarly denote by 4 the nearest neighbor site to 3 
with a connecting bond J43. If J23 is much larger then 
the connecting bonds the spin pair (2, 3) acts as an effec- 
tive singlet. It follows that the strongly correlated singlet 
pair can be frozen out. Due to the virtual triplet excia- 
tions, an effective coupling J14 is generated between the 
sites 1 and 4. These two sites become nearest neighbors 
once the singlet has been eliminated. In a second-order 
perturbation calculation one obtains 
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The new coupling is thus smaller than any of the original 
ones. The energy scale is continuously reduced upon 
iterating the procedure and at the same time the prob- 
ability distribution of the couplings P(J, f2) approaches 
a limiting function. In a gapless random system f2 tends 
to zero at the fixed point of the transformation and the 
low energy tail of the distribution is typically given by: 
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The dynamical exponent z determines how the length 
scale L scales with the time scale r: 



(3.3) 



In general z is not a universal quantity: its value de- 
pends on the form of the original disorder. However z 
stays invariant under renormalizatior£3. Therefore one 
can de duc e its value from the renormalized distribution 
in Eq. ( |3.2| ) . Varying the parameters of the initial distri- 
bution one can reach a situation where the width of the 



distribution in Eq.(3.2) grows without limits, i.e. z for- 
mally tends to infinity. In this case, according to exact 
results on the random AF spin-1/2 chaino, one should 
formally replace z in Eq.( |3.2| ) by — InO, so that the scal- 
ing relation in Eq.(3.3) takes the form: 
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This type of fixed point, where the ratio of any two neigh- 
boring bonds typically tends to zero or infinity, is called 
an infinite randomness fixed point (IRFP). It has been 
conjecture d tha t the MDH renormalization group trans- 
formation (3.1) leads to exact results regarding the sin- 
gular proper ties of the transformation, namely the value 
of tp in Eq. (3.4) is exactE 
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FIG. 4. Schematic RG phase diagram of the random 
dimerized s=l/2 chain as a function of the quantum con- 
trol parameter 8 dim and the energy scale Q. Along the RG 
trajectories the dynamical exponent z(5dim) is constant. The 
non-dimerized model with &dim — is attracted by the IRFP 
with 1/z = 0. With decreasing the energy scale Q, disorder 
in the system is increasing. 
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For the_candom AF spin- 1/2 chain, according to ex- 
act resultsEj any amount of disorder is sufficient to drive 
the system into the IRFP. Similarly, the gapped dimer 
phase will turn into a gapless random dimer phase for 
any amount of disorder, where the dimerization parame- 
ter is defined as 
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in terms of the couplings J Q dd and J even at odd and 
even sites, respectively. The random dimer phase is a 
quantum version of the Griffiths-phase, which has beep 
originally introduced for classical disordered systemsEa. 
The schematic RG-flow diagram of the random dimer- 
ized chain is drawn in Fig. |4]. 




IRFP 

FIG. 5. Schematic RG phase diagram of the random AF 
spin-1 chain, as a function of the disorder strength D and the 
energy scale Q. For weak disorder, D < Do, there is a Hal- 
dane ("H") gap in the spectrum. For intermediated disorder, 
Do < D < Di the system is in the gapless Haldane ("GH") 
phase with a varying dynamical exponent z(D). For strong 
enough disorder, D > Di the system is in the random singlet 
phase and scales into the IRFP. 

The MDH method has also been used to study the sin- 
gular properties of the random AF spin-1 chain. Here 
we first note that the spectrum of the pure system has a 
Haldane gap, which is stable against weak randomness. 
Consequently the MDH renormalization, which is by def- 
inition a strong disorder approach, becomes valid if the 
initial disorder is increased over a limit, say Do- Our 
second remark concerns the "projecting onto the lowest 
level" procedureE3 when after decimating out a strongly 
coupled singlet the generated new coupling is in the form 
of Eq.(3.1), however with a constant of k(S = 1) = 4/3. 
Consequently at the initial RG steps the energy scale 
could behave non-monotonically, so that an IRFP be- 
havior is expected only for strong enough initial disor- 
der. To this problem a modified RG scheme was 
proposedEicS, which is based on the principle of "pro- 



jecting out the highest level". By this method an ef- 
fective Hamiltonian with spin-1 and spin- 1/2 degrees of 
freedom has been introduced, where between the spins 
both AF and F couplings could be present but their 
distribution should respect some constraints. The RG 
analysis of this effective model leads to two different 
types of strongly disordered phases, provided the dis- 
order of the original distribution exceeds the limiting 
value of Do- For an intermediate range of disorder, so 
that D < D < Di the system scales into a quan- 
tum Griffiths-phase, the so-called gapless Haldane phase, 
where z — z(D) is a monotonously increasing function of 
disorder and l/z(D) > 0. When the strength of disorder 
exceeds a second limit, say D > D\ > Do, the dynamical 
exponent becomes infinite and the singular behavior of 
the system is controlled by the IRFP (see Fig ||) . 

Till now there are no numerical estimates about the 
limiting disorder strengths, Dq and D±. For the uniform 



distribution, which corresponds to D = 1 in Eq.(1.2).Jjhc 
system is in the gapless Haldane phase with z « 1.5ej_ 
Finally, we mention the work by Westerberg et al.Eil 
about renormalization of spin-1/2 Heisenberg chains with 
mixed F and AF couplings. In this problem, due to the 
presence of strong F bonds under renormalization spin 
clusters with arbitrary large effective moment S e ^ are 
generated, such that S e ^ ~ A - ", where A is the largest 
local gap in the system and A — > at the fixed point. 
Singularities of different physical quantities are related 
to the scaling exponent u>. 



IV. RENORMALIZATION OF AF SPIN 
LADDERS 

With a ladder geometry, spins are more interconnected 
than in a chain, which leads to a modification of the dec- 
imation procedure described in the previous section. As 
shown in Fig. ^| both spins of a strongly coupled pair, 
say (2,3), are generally connected to the nearest neigh- 
bor spins, denoted by 1 and 4. After decimating out the 
singlet pair the new, effective coupling between 1 and 4 
is of the form: 
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which should replace Eq.(3.1) obtained in the chain 
topology, i.e. with J13 = J42 = 0. With the rule in 
Eq.(4.1) F couplings are also generated. As a conse- 
quence, the renormalized Hamiltonian contains both AF 
and F bonds. When at some step of the renormalization 
an F bond becomes the strongest one, it will lead to the 
formation of an effective spin-1 cluster. In further RG 
steps the system renormalizes into a set of effective spin 
clusters having different moments and connected by both 
AF and F bonds. The detailed renormalization rules have 
already been given in RefO. 
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iteration calc ulated with the disorder parameter D = 1 
(see Eq.(|l.2|)). As seen in the figure the small energy 
tail of the distribution follows the functional form given 
by Eq. (3.2) and the dynamical exponent z given by the 
asymptotic slope of the distributions is finite and has only 
a very weak size dependence. 



FIG. 6. 

geometry. 



Singlet formation and decimation in the ladder 



Conventional ladder, D=1 



chains of the RTI 
for the 2D RFP^ 2 




Due to the ladder topology and the complicated renor- 
malization rules the RG equations can not be treated 
analytically and one resorts to numerical implementa- 
tions of the renormalization procedure. We note that 
a variant of the MDH renormalization has been success- 
fully applied numerically for the twoTdiaacnsional random 
transverse-field Isin#,model (RTIM)c3c3 (also for double 
. An IRFP has been obtained-both 
and the double chain RFIME3. 
In practice we use a finite size version of the MDH 
renormalization, as for the RTIM in Ref [h| In this 
method we start with a finite ladder of L sites and with 
periodic boundary conditions and perform the decima- 
tion procedure until one spin pair with a first gap A re- 
mains in the system. Since A plays the role of the energy 
scale at length scale L, A and L should be related by the 
relation involving the dynamical exponent z. Per- 
forming the above decimation for different samples the 
probability distribution of A in the small A limit is de- 
scribed by the form in Eq.(3.2), where the energy scale 
ft is replaced by L~ z . 

The IRFP is signalled by a diverging z, or more pre- 
cisely the Pt(A)dA distributions have strong L depen- 
dence, so that the appropriate scaling combination is 
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which can be obtained from Eq.(3.2) by formally setting 
z~-lnA~L^. 

In the actual calculations we have considered several 
100.000 realizations of random ladders of length up to 
L = 512. Then, from the distribution of the gap at the 
last step of the RG iteration we have calculated the dy- 
namical exponent, z. The random coup ling s were taken 
from the power-law distribution in Eq.( |l.2| ), where the 
strength of disorder is measured by the parameter D. In 
the following we present our results for the specific ladder 
models discussed in Sec. 2. 



A. Random conventional ladders 

We start with the conventional ladders in Fig. [|.a 
where the couplings along the chains (Jf, r = 1,2) and 
the couplings along the rungs (J/*) are taken from the 
same random distributions. In Fig. ^ we show the prob- 
ability distribution of the gaps at the last step of the RG 
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FIG. 7. Probability distribution of the first gaps for the 
con vent ional random ladder with a disorder D = 1 (see 
and system size L = 32, L = 64, and L = 256. 
For clarity, we have not shown the data corresponding to 
L = 128. The solid lines represent the best fit to the form 
log 10 [P[-log 10 A]] = Al- ^log 10 A, with A 32 = 3.18, 
A e>4 = 2.30, A 12S = 1.92, Aase = 2.14 and z 32 = 0.65, 
264 = 0.90, Zi28 = 1.06, and 2256 = 1.07. We deduce that the 
asymptotic value of the dynamical exponent is ~ 1.07. 



Repeating the calculation for other values of D we ob- 
tain a set of D-dependent dynamical exponents which 
are represented on Fig. |[ For strong disorder we ob- 
tain z(D) < D, which means that disorder is reduced 
in the course of the-renormalization. In the terminology 
of Motrunich et aiJlj, this system is a Unite randomness 
system, as opposed to the infinite randomness systems 
that will be considered in the next subsections. For weak 
disorder the dynamical exponent predicted by the ap- 
proximative MDH renormalization is lowered below one 
for D < Do » 1. Here we argue that in this region 
the effect of disorder is irrelevant, so that the system is 
in the gapped RUS phase. Indeed, in a pure quantum 
system, where scaling in time and space is isotropic, the 



dynamical exponent is z pure = 1. Similarly, for disor- 
der induced gapless systems, where disorder in the time 
direction is strictly correlated, the dynamical exponent 
can not be smaller, than in the pure system, so that 

Zdis ^ Zp Ur 



= 1. Consequently, if the disorder induced 
dynamical exponent is z^ < 1, then disorder could only 
influence the correction to scaling behavior, but the sys- 
tem stays gapped. In view of this remark Dq can be 
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considered as the lower limiting value of the disorder, 
where the conventional finite randomness behavior ends. 
So that the phase diagram of random conventional two- 
leg spin ladders consists of two phases: a gapped RUS 
(Haldane) phase and a random gapless Haldane phase. 
The latter is characterized by a finite dynamical expo- 
nent z(D) for any strong but finite initial disorder. Con- 
sequently there is an important difference with the ran- 
dom AF spin-1 chain which flows into the IRFP above a 
finite critical value of randomness (see Fig. |J) • 
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FIG. 8. Variation of dynamical exponent z versus disorder 
D for the conventional ladder with lengths L = 32, L = 64, 
L = 128, and L — 256. For large system sizes and strong 
disorder, one has Zoa ~ 0.42 + 0.58D < D. In the region with 
Zoo < 1) where the disorder is irrelevant, the system is in the 
gapped Haldane phase (see text). 



far from the maximum the convergence of the data is 
faster. To decide about the possible limiting value of z, 
in particular at the maximum of the curves, we analyze 
the behavior for 7 = 1 in Fig. |^c, which is just a dimer- 
ized random chain, the .properties of which are exactly 
known by some extentc3E3. 
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FIG. 9. Finite-size estimates of the dynamical exponent of 
random conventional ladders with staggered dimerization as 
a function of the coupling ratio, Jmax/ Jmax, with a disorder 
parameter D = 1 and for different dimerizations: a) 7 = 0.5, 
b) 7 = 0.75, c) 7 = 1.. In d) a similar calculation for the 
random XX-chain is presented (D — 1, 7 = 1), where the 
exact dynamical exponent in Eq.(4.3), obtained in the L — > 00 
limit, is shown by dashed line. 



B. Random ladders with staggered dimerization 

In this subsection we consider conventional ladders 
with staggered dimerization having a dimerization pa- 
rameter, < 7 < 1, in Eq.(2.2). The different type 
of couplings in the ladder are taken from the power-law 
distribution in Eq.( |1.2| ), each having the same disorder 
parameter, D, however the range of the distribution for 
the different type of couplings are, < Jf 1 < J^ ax 
for the rung couplings, < J™ eak < (1 — 7) J max and 
< j ; itro ™ 9 < (1 + j)J max for the weaker and stronger 
chain couplings, respectively. For a fixed value of 7 and 
D we have calculated the finite-size dependent effective 
dynamical exponent, function of the coupling ra,- 

J max I Jmax- 

As shown in Fig. || the effective exponents have the 
same type of qualitative behavior for different values of 
the dimerization parameter, 7. In each case the curves 
have a maximum at some value of the couplings, where 
the finite-size dependence is the strongest, whereas more 



The random critical point of this system is situated 
a * J max I ^rnax = 2, where the critical behavior is gov- 
erned by an IRFP, so that the dynamical exponent, z, is 
formally infinity. For any other values of the couplings, 
Jmax/ Jmax, the system is in the random dimer phase, 
where the dynamical exponent is finite and coupling de- 
pendent. To see the general tendency of finite-size con- 
vergence of the z exponent around the critical point we 
have repeated the calculation at the D = 1, 7 = 1 case 
for the random AX-chain and the numerical finite-size 
results are compared in Fig. ^d with the exact value of 
the dynamical exponent, as given by the solution of the 
equation: 



J R 



J rt 



D 2 



D 2 



(4.3) 



known from Ref. |3^j3^. As seen in Fig. ^| the dynam- 
ical exponents of random XX- and Heisenberg-ladders 
have very similar coupling dependence and one expects 
the same type of divergence at the critical point for all 
values of 7. In Fig. 10 we illustrate the scaling behavior 
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of the gap at the transition point, i.e. at the maximum 
of the curves in Fig. ^|a. 





randomness behavior. For different disorder parameter, 
-D, the position of the random critical line is modified, 
generally stronger disorder is in favor of the random rung 
singlet phase, see Fig. O . 

We note that the previously studied random conven- 
tional ladder is contained as a special point in this phase 
diagram at Jmax/Jmax = 1 and 7 = 0. This point is in 
the random rung singlet phase for any value of D, thus 
the dynamical exponent is finite in accordance with the 
previous results. 



(a) 



(b) 



FIG. 10. a) Probability distribution of the first gap at 
the transition point of the random conventional ladder with 
staggered dimerization, D = 1, 7 = 0.5, Jmax/Jmax = 1.1. 
The distributions become broader and broader with L, which 
signals infinite randomness behavi or, b) Scaling plot in terms 
of the scaling combination in Eq.(3.4). 



oD = 1 
■ D = 0.75 
v D = 0.5 
a D = 0.25 




FIG. 11. Phase diagram of random conventional ladders 
with staggered dimerization for different disorder parameters. 
The random dimer (RD) phase and the random rung singlet 
(RRS) phase are separated by a random critical line of infinite 
randomness behavior. 

The distributions in Fig. [To| a become broder and 
broder with the size and the effective dynamical exponent 
is increasing with the size without limits. An appropri- 
ate scaling collaps of the gap-distributions has bee n ob - 
tained in Fig. [lOb, where the scaling variable in Eq.(|4.2|), 
with ip = 1/2 is used. Similar type of infinite random- 
ness behavior is observed at other points of the critical 
lines, with the same exponent ip — 1/2, which turned 
out to be universal. Thus we conclude that the random 
conventional ladder with staggered dimerization has two 
Griffiths-type gapless phases, the random dimer phase 
and the random rung singlet phase, which are separated 
by a random critical line, along which there is infinite 



C. Random zig-zag ladders 

For the zig-zag ladders the nearest neighbor couplings 
(jf 1 = Jf = Ji) are taken from the power-law distribu- 
tion in Eq. flLj ) with the coupling J\ within the range 
< J\ < Jjn ax - Similarly, the next-nearest neighbor 
couplings (J; = Jf) are taken from the same type of 
power-law distribution and the range of couplings is now 

J 1 < J max- 



zig-zag ladder, D=1 




zig-zag ladder, D=5 

~~ ' ' ' L=128 

L=256 — 
L=512 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.1 

Jmax/Jmax 

(a) 



0.2 0.4 0.6 0.8 1 1.2 

Jmax/Jmax 

(b) 



FIG. 12. Variation of the dynamical exponent z versus 
Jmax/Jmax f° r the zig-zag ladder with a disorder D = 1 (a) 
and D — 5 (b). For D = 1, the MDH renormalization group 
data with L — 128, L — 256 and L — 512 have been com- 
pared to the DMRG calculation with L — 32 (see Fig. p^j ). 
The lines connecting the calculated points are guides to the 
eyes. 

The calculated dynamical exponent, z, as shown in 
Figs. [l2] has its maximum at Jmax/Jmax = an d around 
this point one can observe strong finite-size dependence, 
the range of which is wide, in particular for weak disor- 
der (see Fig. H). At J max /J max = 0, where the zig-zag- 
ladder reduces to a random AF chain, the system is in 
the IRFP, thus the extrapolated value of the dynamical 
exponent is formally infinity. Given the strong finite-size 
corrections in the numerical RG data of the dynamical 
exponentcj in Fig. 12. it is difficult to decide whether the 
IRFP behavior of the zig-zag ladders is extended to a fi- 
nite region of the couplings J max /J max > or whether 
this region shrinks to a single point only. The first sce- 
nario may be related to the existence of a gapless phase 
of the pure model for J 2 / J 1 < 0.24. 
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FIG. 13. Distribution of the log surface magnetization 
of the random zig-zag ladder using a power-law distribution 
(D = 1, Jmax/Jmax = 0.5), for different lengths of the ladder, 
L. The asymptotic slope of the distribution, indicated by the 
stra ight line, is the inverse of the dynamical exponent, see 
Eq.Q. 



To discuss this issue we have calculated the dynamical 
exponent by an independent method based on density 
matrix renormalization. In principle, the dynamical ex- 
ponent is related to the distribution of the first gap, A, 
in the small A limit, see in Eq.(3.2) with A — * J. How- 
ever, a precise numerical calculation of a small A by the 
DMRG method is very difficult, therefore we used an- 
other strategy, as described in details in Refs. H?],|d| By 
this method one considers the equivalent AF chain with 
random first- and second-neighbor couplings (see Fig. 
Id) and with fixed-free boundary conditions and calcu- 
late the surface magnetization, m s , at the free end, which 
can be done very accurately by the DMRG method. As 
argued in Refs. [5^,^4] for a random chain m s and A can 
be considered as dual quantities, so that the distribution 
of the surface magnetization is asymptotically given by 



P(hxm s ) 



0. 



(4.4) 



Thus the dynamical exponent z can be obtained from 
an analysis of the small m s tail of the distribution, as 
illustrated in Fig. [l3| where the distribution function of 
lnm s is given in a log-log plot for different lengths of 
the ladder. As seen in this Figure the slope of the dis- 
tribution is well defined for larger systems, from which 
one can obtain an accurate estimate for the dynamical 
exponent, which is finite. Repeating the calculation for 
other values of the coupling ratio, J max / J max i we have- 
obtained a set of the dynamical exponents, which are 
plotted in Fig. [l2| These accurate DMRG data show 
that the extrapolated values of the effective exponents 
calculated by the numerical RG- method arc finite for any 
J max /J max > 0- Consequently the random zig-zag lad- 
der has just one IRFP at J max I J max = 0; whereas the 



system in the region of J max /J max > is in a gapless 
random dimer phase. In view of the numerical results 
in Fig. [l^, where z,^ seems to stay over z pure = 1, it is 
quite probable that the randum dimer phase exists for 
any small value of the disorder. 



D. Random J\—Ji ladders 



The full ladder, as represented in Fig. [l].e has three 
different type of couplings: Ji, jf- and Jf. Here we con- 
sider a special case of this model, when the chain (Ji) and 
rung (Jp) couplings are taken from the same power-law 
distribution with a disorder parameter D and having a 
range of < <// , Jf 1 < J max . On the other hand the diag- 
onal couplings are taken from the same type of power-law 
distribution and are within the interval < Jf < J max - 
This model, having first- and second-neighbor interac- 
tions, is called a J\-Ji ladder. We have calculated the 
finite-size dynamical exponents as a function of the cou- 
pling ratio Jmax/Jmax f° r different strengths of disorder 
(see Fig. ||). These curves show similar qualitative be- 
havior as those calculated for the random conventional 
ladders with staggered dimerization in Figs. ||, so that 
we can draw similar conclusions. 



J1-J2 ladder, D=1 



J1-J2 ladder, D=5 




| 4 

1 3 
O 2 

I 






L=32 — i — 
L=64 — x— 
L=128 - 
L-256 a 




H. 
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,„ -5K D 1 











0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 
(a) 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

(b) 



FIG. 14. Variation of the dynamical exponent z versus 
Jmax/Jmax f° r the random J\-Ji ladder with a disorder 
D — 1 (a) and D — 5 (b). The lines connecting the calcu- 
lated points are guides to the eyeSf-Note the strong finite-size 
corrections in the Griffiths phasesCj. 

The extrapolated position of the maximum of the z 
curves is identified as a quantum critical point with in- 
finite randomness behavior. Indeed, repeating the cal- 
culation as indicated for the dimerized ladder mod el i n 
Fig. |h] we obtained a scaling behavior as in Eq.(3.4). 
with an exponent which is compatible with ip = 1/2. 
The random quantum critical point separates two gapless 
Haldane phases, having odd and even topological order, 
respectively. 
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0.6 



0.5 



if 0.4 



0.3 - 




0.2 



0.1 



Disorder D 

FIG. 15. Phase diagram of the random J\—J% ladder ob- 
tained from the maximum in the variation of the dynami- 
cal exponent versus J max I J max for the largest available sizes 
L = 256 (see Fig. ^ for D = 1 and D = 5). The two straight 
lines connecting the calculated transition points are guides to 
the eyes. The whole transition line between the two phases 
with even and odd topological orders, respectively, is presum- 
ably a line of IRFPs. 

Repeating the calculation for different disorder param- 
eters we obtain a phase diagram shown in Fig. [l5| In 
the range of disorder we used in the calculation the ran- 
dom critical point is always attracted by the IRFP, this 
property probably remains true for any small value of 
disorder. 



V. DISCUSSION 

In this paper different type of random AF spin ladder 
models have been studied by a numerical strong disor- 
der RG method. In particular we asked the questions i) 
how the phase diagrams of the pure models are modified 
due to quenched disorder and ii) how the concepts ob- 
served in random AF chains, such as infinite randomness 
and Griffiths-type singularities are valid for these more 
complicated, quasi-one-dimensional models. 

In our numerical calculations we observed as a general 
rule that for strong enough disorder the ladder models, 
like the random chains, become gapless. The dynamical 
exponent of the models is generally non-universal: z de- 
pends on both the strength of disorder and on the value 
of the couplings. In models where there is a competition 
between different types of phases, either due to staggered 
dimerization or due to frustration, such as in the J1-J2 
model, at the phase boundary the critical behavior of 
the random model is generally controlled by an infinite 
randomness fixed point, at least for strong enough dis- 
order. The low-energy properties of the systems in this 



IRFP are asymptotically exactly knowrufrjarn analytical 
calculations in random AF spin chainsu'El Thus the 
general phase diagram consists of Griffiths-type phases 
with different topological order separated by a random 
critical point of the IRFP type. The zig-zag ladder is an 
exception, where there is just one Griffiths phase and the 
random critical point is located at its boundary. 

Next we turn to discuss about possible cross-over ef- 
fects when the strength of disorder is varied. These prob- 
lems can not be directly studied by the simple strong dis- 
order RG method, however, from arguments considering 
the sign of z^is — z pure and from analogous investigations 
on quantum spin chainsEiHa we can suggest the follow- 
ing picture. Originally gapped phases could stay gapped 
for weak disorder and become gapless only if the strength 
of disorder exceed some finite limiting value, as seen for 
the random conventional ladder. However, for frustrated 
ladders, such as the zig-zag and the J\ — J 2 ladders, any 
small amount of disorder seems to bring the system into 
a random gapless phase. At a phase boundary, such as in 
the staggered dimerized ladder and the J\—Ji model, the 
random critical behavior is of the IRFP type, probably 
for any small amount of disorder. 

At this point we comment on the similarity of the low- 
energy behavior of spin chains with 2S = odd (2S = 
even) spins and that of spin-1/2 ladders with n = odd 
(n = even) legs. If the pure systems are gapless, i.e. 
2S = n = odd, strong enough disorder is expected to 
bring both systems into the IRFP. For 2S = n = odd > 3, 
there is a limiting disorder strength, D c (n), below which 
the system is described by a conventional random fixed 
point with z < 00. On the other hand for 25* = n = even 
we have only a partial analogy: for weak disorder both 
systems are gapped, which turn into a gapless Griffiths- 
type phase for stroger disorder. While the random ladder 
stays in this conventional random phase for any strength 
of the disorder the random spin chain will turn into a 
IRFP behavior at some finite limiting randomness. This 
type of infinite randomnass behavior can, however, be 
seen for frustrated even-leg J\ — Ji ladders at the tran- 
sition point. We can thus conclude that random ladders 
with even and odd number of legs belong to different 
universality classes. 

Finally, we comment on random square lattice anti- 
ferromagnets, which can be obtained in the limit when 
the number of legs, n, goes to infinity. By increasing n 
the value of the limiting strength, D c (n), is expected to 
increase, too, both for n = even (for frustrated ladders) 
and n = odd. In the limit, n — > 00, D c (n) very probably 
tends to infinity, so that the critical behavior of that sys- 
tem is described by a conventional random fixed point. 
Work is in progress to verify this scenario and to obtain a 
general physical picture about the low-energv-properties 
of random two-dimensional antiferromagnetsc3. 
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